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Abstract 



Charged-hadron transverse-momentum and pseudorapidity distributions in proton- 
proton collisions at v^ = 7 TeV are measured with the inner tracking system 
of the CMS detector at the LHC. The charged-hadron yield is obtained by count- 
ing the number of reconstructed hits, hit-pairs, and fully reconstructed charged- 
particle tracks. The combination of the three methods gives a charged-particle multi- 
plicity per unit of pseudorapidity dNjh/'^'/||n|<o.5 = ^-^^ ^ 0-01 (stat.) ± 0.23 (syst.) 
for non-single-diffractive events, higher than predicted by commonly used mod- 
els. The relative increase in charged-particle multiplicity from ^s = 0.9 to 7 TeV 
is (66.1 ± 1.0 (stat.) ± 4.2 (syst.))%. The mean transverse momentum is measured to 
be 0.545 ± 0.005 (stat.) ± 0.015 (syst.) GeV/c. The results are compared with similar 
measurements at lower energies. 
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Introduction. Measurements of particle yields and kinematic distributions are an essential first 
step in exploring a new energy regime of particle collisions. Such studies contribute to our 
understanding of the physics of hadron production, including the relative roles of soft and 
hard scattering contributions, and help construct a solid foundation for other investigations. 
In the complicated environment of LHC pp collisions HI, firm knowledge of the rates and 
distributions of inclusive particle production is needed to distinguish rare signal events from 
the much larger backgrounds of soft hadronic interactions. They will also serve as points of 
reference for the measurement of nuclear-medium effects in Pb-Pb collisions in the LHC heavy 
ion program. 

The bulk of the particles produced in pp collisions arise from soft interactions, w^hich are mod- 
eled only phenomenologically Experimental results provide the critical guidance for tuning 
these widely-used models and event generators. Soft collisions are commonly classified as 
elastic scattering, inelastic single-diffractive (SD) dissociation, double-diffractive (DD) dissoci- 
ation, and inelastic non-diffractive (ND) scattering ^. (Double-Pomeron exchange is treated 
as DD in this paper.) All results presented here refer to inelastic non-single-diffractive (NSD) 
interactions, and are based on an event selection that retains a large fraction of the ND and DD 
events, while disfavouring SD events. 

The measurements focus on transverse-momentum pj and pseudorapidity rj distributions. The 
pseudorapidity commonly used to characterize the direction of particle emission, is defined as 
f] = — lntan(0/2), where 9 is the polar angle of the direction of the particle with respect to 
the anti-clockwise beam direction. The count of primary charged hadrons, N^-hz is defined to 
include decay products of particles with proper lifetimes less than 1 cm. Products of secondary 
interactions are excluded, and a per-cent-level correction is applied for prompt leptons. The 
measurements reported here are of dNch/drj and dN^h/dpr in the pseudorapidity range \rj\ < 
2.4 and closely follow our previous analysis of minimum-bias data at lower centre-of-mass 
energies of a/s = 0.9 and 2.36 TeV as reported in Ref. |3|. 

The data for this study are drawn from an integrated luminosity of 1.1 f/b^^ recorded with the 
Compact Muon Solenoid (CMS) experiment [4J on March 30th, 2010, during the first hour of 
the LHC operation at^/s = 7 TeV. These results are the highest centre-of-mass energy measure- 
ments of the dN^\^/dr] and dN^^\^/dpj distributions conducted at a particle collider so far and 
complement the other recent measurements of the ALICE experiment at 7 TeV |i5J. 

Experimental methods. A detailed description of the CMS experiment can be found in Ref. f^. 
The detectors used for the present analysis are the pixel and silicon-strip tracker (SST), covering 
the region \f]\ < 2.5 and immersed in a 3.8 T axial magnetic field. The pixel tracker consists 
of three barrel layers and tw^o end-cap disks at each barrel end. The forward calorimeter (HE), 
which covers the region 2.9 < \t]\ < 5.2, was also used for event selection. The detailed Monte 
Carlo simulation (MC) of the CMS detector response is based on Geant4 [6]. 

The event selection and analysis methods in this paper are identical to those used in Ref. (Sj, 
where more details can be found. The inelastic pp collision rate was about 50 Hz. At these 
rates, the fraction of events in the data, where two or more minimum-bias collisions occurred 
in the same bunch crossing, is estimated to be less than 0.3% and w^as neglected. Any hit in the 
beam scintillator counters (BSC, 3.23 < \f]\ < 4.65) coinciding with colliding proton bunches 
was used for triggering the data acquisition. A sample mostly populated with NSD events 
was selected by requiring a primary vertex (PV) to be reconstructed with the tracker, together 
•with at least one HE tower in each end with more than 3 GeV total energy. Beam-halo and 
other beam-background events were rejected as described in Ref. [3J. The remaining fraction 
of background events in the data w^as found to be less than 2 x 10^^. The numbers of events 



Table 1: Numbers of events passing the selection cuts. The selection criteria are applied in 
sequence, i.e., each line includes the selection from the previous ones. 



Selection Number of events 
Colliding bunches + one BSC signal 68 512 

Reconstructed PV 61 551 

HF coincidence 55 113 

Beam-halo rejection 55 104 

Other beam-background rejection 55 100 



Table 2: Fractions of SD, DD, ND, and NSD processes obtained from the PYTHIA and PHO- 
JET event generators before any selection, and the corresponding selection efficiencies deter- 
mined from the MC simulation. 



PYTHIA 


PHOJET 




Frac. 


Sel. eff. 


Frac. Sel. eff. 


SD 


19.2% 


26.7% 


13.8% 30.7% 


DD 


12.9% 


33.6% 


6.6% 48.3% 


ND 


67.9% 


96.4% 


79.6% 97.1% 


NSD 


80.8% 


86.3% 


86.2% 93.4% 



satisfying the selection criteria are listed in Table [T] 

The event selection efficiency was estimated with simulated events using the PYTHIA [Tl |8l 
and PHOJET ll9l [T0ll event generators. The relative event fractions of SD, DD, and ND processes 
and their respective event selection efficiencies are listed in Table |2] The fraction of diffractive 
events is predicted by the models to decrease as a function of collision energy, while the selec- 
tion efficiency increases. At a/s = 7 TeV, the fraction of SD (DD) events in the selected data 
sample, estimated with PYTHIA and PHOJET, are 6.8% (5.8%) and 5.0% (3.8%), respectively, 
somewhat higher than at a/s = 0.9 and 2.36 TeV HS). With PYTHIA, the overall correction 
for the selection efficiency of NSD processes and for the fraction of SD events remaining in 
the data sample lowers the measured charged-particle multiplicity by 6% compared with the 
uncorrected distribution. 

The dNf.]^/dri distributions were obtained, as in Ref. [3|, with three methods, based on count- 
ing the following quantities: (i) reconstructed clusters in the barrel part of the pixel detector; 
(ii) pixel tracklets composed of pairs of clusters in different pixel barrel layers; and (iii) tracks 
reconstructed in the full tracker volume. The third method also allow^s a measurement of the 
dNch/dpi distribution. All three methods rely on the reconstruction of a PV [11]. The PV re- 
construction efficiency was found to be 98.3% (98.0%) in data (MC), evaluated after all other 
event selection cuts. In case of multiple PV candidates, the vertex w^ith the largest track mul- 
tiplicity w^as chosen. The three methods are sensitive to the measurement of particles dow^n to 
Pt values of about 30, 50, and 100 MeV/c, respectively. Only 0.5, 1.5, and 5% of all charged 
particles are estimated to be produced below these pj values, respectively, and these fractions 
were corrected for. 

The measurements were corrected for the geometrical acceptance (~2%), efficiency (~5-10%), 
fake (<1%) and duplicate tracks (<0.5%), low-px particles curling in the axial magnetic field 
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Figure 1: Differential yield of charged hadrons in the range \t]\ < 2.4 in 0.2-unit-wide bins of 
I?/ 1 in NSD events. The solid curves represent fits of Eq. n^to the data. The measurements with 
increasing rj are successively shifted by six units along the vertical axis. 

(<1%), decay products of long-lived hadrons (<2%) and photon conversions (<1%) and in- 
elastic hadronic interactions in the detector material (?»l-2%), where the size of the corrections 
in parentheses refer to the tracking method. The PYTHIA parameter set from Ref. |8| was 
chosen to determine the corrections, because it reproduces the dN^-^^/di] and charged-particle 
multiplicity distributions, as well as other control distributions at 7 TeV, better than other avail- 
able tuning parameter sets. Although the corrections do not depend significantly on the model 
used, it is indeed important that the simulated data set contains a sufficient number of high- 
multiplicity events to determine these corrections with the desired accuracy. 

Results. For the measurement of the dN^\y/ dipj distribution, charged-particle tracks with pj in 
excess of 0.1 GeV/c were used in 12 different | j/| bins, from to 2.4. The average charged-hadron 
yields in NSD events are shown in Fig. n^as a function of px and \rj\. The Tsallis parametrization 
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was fitted to the data. The pj spectrum of charged hadrons, 1/ {2npj)d^N^h/dri dpj, measured 
in the range \f]\ < 2.4, is shown in Fig. |2] for data at 0.9, 2.36, and 7 TeV. The high-px reach of 
the data is limited by the increase of systematic uncertainties with pj. The fit to the data (Eq.[l]l 
is mainly used for extrapolations to px = 0, but is not expected to give a good description of 
the data in all tj bins with only two parameters. The parameter T and the exponent n w^ere 
found to be T = 0.145 ± 0.005 (syst.) GeV and n = 6.6 ± 0.2 (syst.). The average px, calculated 
from a combination of the measured data points and the low- and high-px contributions as 
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Figure 2: Charged-hadron yield in the range \f]\ < 2.4 in NSD events as a function of pj; the 
systematic uncertainties are smaller than the symbols. The measurements at ^/s = 0.9 and 
2.36 TeV [3J are also shown. The solid lines represent fits of Eq.lTlto the data. 

determined from the fit, is (pj) = 0.545 ± 0.005 (stat.) ± 0.015 (syst.) GeV/c. 

Experimental uncertainties related to the trigger and event selection are common to all the 
analysis methods. The uncertainty related to the presence of SD (DD) events in the final sample 
was estimated to be 1.4% (1.1%), based on consistency checks between data and simulation for 
diffractive event candidates. The total event selection uncertainty, which also includes the 
selection efficiency of the BSC and HF, w^as found to be 3.5%. Based on studies similar to those 
presented in Ref. [3], additional 3% and 2% uncertainties were assigned to the tracklet and track 
reconstruction algorithm efficiencies, respectively. Corrections at the percent level were applied 
to the final results to extrapolate to pj = 0. The uncertainty on these extrapolation corrections 
was found to be less than 1%. All other uncertainties are identical to those listed in Ref. ||3j. The 
dNch/df] measurements were repeated on a separate data sample without any magnetic field, 
for which almost no pj extrapolation is needed, and gave results consistent within 1.5%. The 
final systematic uncertainties for the pixel counting, tracklet, and track methods w^ere found to 
be 5.7%, 4.6%, and 4.3%, respectively, and are strongly correlated. 

For the dNch/drj measurements, the results for the three individual layers within the cluster- 
counting method w^ere found to be consistent w^ithin 1.2% and were combined. The three layer- 
pairs in the pixel-tracklet method provided results that agreed within 0.6% and were also com- 
bined. Finally, the results from the three different measurement methods, which agree w^ith 
the combined result w^ithin 1% to 4% depending on rj, were averaged. The final dNt^^/drj dis- 
tributions are shown in Fig. p^for a/s = 0.9, 2.36, and 7 TeV. The CMS results are compared 
w^ith measurements made by other experiments. In the ATLAS Collaboration analysis [17J, 
events and particles were selected in a different region of phase space, which makes a direct 
comparison difficult. Their results are therefore not included in the figure. 

The results can also be compared to earlier experiments as a function of v^. The energy depen- 
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Figure 3: Distributions of dNch/drj, averaged over the three measurement methods and com- 
pared with data from UA5 [15J (pp, with statistical errors only) and ALICE I.16J (w^ith systematic 
uncertainties). The shaded band shows systematic uncertainties of the CMS data. The CMS and 
UA5 data are averaged over negative and positive values of ?/. 

dence of the average charged-hadron pj can be described by a quadratic function of Ins llTSl . 
As shown in Fig.|4| the present measurement follows this trend. The choice of the \?]\ interval 
can influence the average pj value by a few per cent. 

For \ri\ < 0.5, the average charged multiplicity density is dNf-h/drj = 5.78 ± 0.01 (stat.) ± 
0.23 (syst.) for NSD events. The a/s dependence of the measured dNch/rf?/|j;«o is shown in 
Fig. [5} which includes data from various other experiments. The dN^^^f drj results reported here 
show a rather steep increase betw^een 0.9 and 7 TeV, which is measured to be (66.1 ± 1.0 (stat.) ± 
4.2(syst.))%. Using a somewhat different event selection, the ALICE collaboration has found 
a similar increase of (57.6 ± 0.4 (stat.)^^ g (syst.))% [5J. The measured charged-particle multi- 
plicity is accurate enough to distinguish among most sets of event-generator tuning parameter 
values and various models. The measured value at 7 TeV significantly exceeds the predic- 
tion of 4.57 from PHOJET |9, 10 1, and the predictions of 3.99, 4.18, and 4.34 from the DW Il22l , 
ProQ20 [23 1, and PerugiaO [24J tuning parameter values of PYTHIA, respectively, while it is 
closer to the prediction of 5.48 from the PYTHIA parameter set from Ref. |8| and to the re- 
cent model predictions of 5.58 and 5.78 from Refs. Il25ll and ||26|. The measured excess of the 
number of charged hadrons with respect to the event generators is independent of t] and con- 
centrated in the pj <\ GeV/c range. These differences indicate the need for a continued model 
development and simulation tuning. Work on updated event generators based on LHC data is 
currently underway. 

Summary. Charged-hadron transverse-momentum and pseudorapidity distributions have been 
measured in proton-proton collisions at yjs = 7 TeV. The numerical values of the data pre- 
sented in this paper can be found in the HEPDATA database [|32j . The combined result for 
the central pseudorapidity density, from three mutually consistent methods of measurement, 
is dN^^/drj \ |„|<o.5 = 5. 78 ± 0.01 (stat.) ± 0.23 (syst.) for non-single-diffractive events. This value 
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Figure 4: Average px of charged hadrons as a function of the centre-of-mass energy. The CMS 
measurements are for \r]\ < 2.4. Also shown are measurements from the ISR [ 19 1 {^X')' E735 [20 1 
(pp), and CDF im (pp) for |?/ 1 < 0.5, and from UAl lUlI {ff) iox\y]\ < 2.5. The solid line is a fit 
of the functional form (px) = 0.413 - 0.0171 In s + 0.00143 In^ 
the CMS data include the systematic uncertainties. 



s to the data. The error bars on 



is higher than most predictions and provides new information to constrain ongoing improve- 
ments of soft particle production models and event generators. The mean transverse momen- 
tum has been measured to be 0.545 ± 0.005 (stat.) ± 0.015 (syst.) GeV/c. These studies are the 
first steps in the exploration of particle production at the new centre-of-mass energy frontier, 
and contribute to the understanding of the dynamics in soft hadronic interactions. 
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